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We investigate the influence of exciton-phonon coupling on the dynamics of a strongly coupled 
quantum dot-photonic crystal cavity system and explore the effects of this interaction on different 
schemes for non-classical light generation. By performing time-resolved measurements, we map 
out the detuning-dependent polariton lifetime and extract the spectrum of the polariton-to-phonon 
coupling with unprecedented precision. Photon-blockade experiments for different pulse-length and 
detuning conditions (supported by quantum optical simulations) reveal that achieving high-fidelity 
photon blockade requires an intricate understanding of the phonons’ influence on the system dy¬ 
namics. Finally, we achieve direct coherent control of the polariton states of a strongly coupled 
system and demonstrate that their efficient coupling to phonons can be exploited for novel concepts 
in high-fidelity single photon generation. 


INTRODUCTION 

The strong coupling between a single photon and a 
single quantum emitter is of substantial interest for both 
investigations of the fundamentals of quantum optics as 
well as potential applications in optical computing, quan¬ 
tum metrology, and quantum cryptography [1-3]. This 
universality is reflected by the diversity of associated 
experimental realizations, ranging from those in atomic 
physics [4-6] to superconducting systems [7, 8] and semi¬ 
conductor devices [9]. In the solid state, self-assembled 
quantum dots (QDs) are the most investigated quantum 
emitters due to their strong interaction with light as well 
as their nearly transform-limited linewidth [10-12]. Opti¬ 
cal resonators, in the form of both micropillar structures 
[9] and photonic crystal cavities [13], are widely used to 
enhance this light-matter interaction. Photonic crystal 
cavities are especially promising for on-chip integration 
of quantum optical circuits due to the convenient fabrica¬ 
tion of integrated waveguide and detector structures [14]. 
In contrast to other systems, semiconductor quantum 
emitters are usually embedded in a crystalline host ma¬ 
trix, resulting in strong interactions with phonons. While 
in many circumstances such interactions are undesirable, 
it is also possible to carefully design experiments which 
benefit from the presence of phonons. For example, it was 
recently demonstrated that the QD coupling to phonons 
can be exploited for a robust and high-fidelity prepara¬ 
tion of exciton and biexciton states [15-17]. Meanwhile, 
when QDs are coupled to semiconductor optical cavities, 
far off-resonant feeding effects have been shown to be 
consistent with a phonon-induced interaction [18, 19]. In 
fact, for QD-cavity systems the coupling to phonons and 


its impact on applications has been extensively studied in 
the weakly coupled regime [20, 21]. However, for strongly 
coupled systems these effects remain largely unexplored 
experimentally. As such systems become increasingly rel¬ 
evant for their potential in building quantum optical net¬ 
works, it is essential to develop a solid understanding of 
the influence that phonons have on the system dynamics. 

In this paper, we investigate the interaction of a 
strongly coupled QD-cavity system with phonons and 
explore the impact of this coupling on different schemes 
for non-classical light generation. By performing tem¬ 
porally resolved spectroscopy experiments, we map out 
the detuning-dependent polariton lifetime and extract 
the spectrum of the polariton-to-phonon coupling. As 
a result, we conclude that the high efficiency of this 
coupling entails its importance for all applications of 
strongly coupled solid-state systems. To corroborate this 
finding, we investigate the influence of electron-phonon 
interaction on one of the key applications of strongly 
coupled systems: non-classical light generation. To 
this end, we perform photon-blockade experiments and 
complementary quantum optical simulations under 
different pulse length and detuning conditions, which 
reveal that the pulse length has to be chosen correctly 
for a given detuning in order to obtain high-fidelity 
photon blockade. In particular, the impact of coupling 
to phonons is found to be more pronounced for the case 
of detuned photon blockade, which is the condition that 
yields the highest fidelity and efficiency [22]. Finally, we 
demonstrate that the efficient coupling to phonons can 
also be utilized for novel concepts in high-fidelity single 
photon generation. 
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STRONGLY COUPLED QD-CAVITY SYSTEMS 


The sample investigated consists of a layer of low den¬ 
sity InAs QDs grown by molecular beam epitaxy and 
embedded in a photonic crystal L3 cavity [23]. The en¬ 
ergy structure of a QD strongly coupled to an optical 
cavity (without phonon coupling) is usually described by 
the Jaynes-Cummings (JC) Hamiltonian 


H = LOaCL 1 a + (u„ + A )<rV + g(a) <7 + acr^) (1) 


where tu a denotes the frequency of the cavity, a the cavity 
mode operator, a the lowering operator of the quantum 
emitter, A the detuning between quantum emitter and 
cavity, and g the emitter-cavity field coupling strength. 
After introducing dissipation into the JC system we ob¬ 
tain a quantum Liouville equation with complex eigenen- 
ergies [24]: 
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FIG. 1. (a) Schematic illustration of the energy structure of 
a strongly coupled system as given by the Jaynes-Cummings 
ladder obtained from equation 2. The blue arrow illustrates 
resonant excitation of UP1, the red arrows photon decay 
and the purple arrow phonon-assisted transfer, (b) Cross- 
polarized reflectivity spectra of the coupled QD-cavity system 
obtained by temperature-tuning the QD through the cavity 
resonance. An anticrossing of the peaks clearly demonstrates 
the strong coherent coupling. 


where E± corresponds to the nth rung of the system and 
k and 7 to the cavity and QD energy decay rates, re¬ 
spectively. The real parts of E± yield the energies of the 
states whereas the imaginary parts yield their linewidths. 
Note that all four parameters k, 7 , A and g contribute 
to the polariton splittings and linewidths. 

The resulting eigenenergies, the Jaynes-Cummings- 
ladder dressed states, are presented in figure la. They 
consist of a series of anticrossing branches that are la¬ 
belled UPn (LP?r) for the upper (lower) polariton where 
n is the index of the rung. For n photons in the cavity the 
energy is noj a (dashed red lines), and the energy of the 
quantum emitter (dotted green line) varies with a detun¬ 
ing parameter. Due to the coupling, the resulting energy 
eigenstates are the anticrossing polariton branches. At 
resonance, the splitting is given by 2 g^/n. 

Experimentally, we can observe the splitting of the 
lowest-manifold polaritons in cross-polarized reflectivity 
[25] through control of the crystal lattice temperature. 
The result of such an experiment is presented in figure lb 
and shows a clear anticrossing. A fit to the data reveals 
values of g = 12.3-27T GHz and k = 18.4-2-7T GHz. The ra¬ 
diative lifetime of a QD in a bulk photonic environment 
is known to be ~ 1 ns and thus, for QDs in photonic 
crystal cavities 7 is much smaller than n. Due to the 
photonic bandgap 7 is further suppressed and reported 
values for the continuum mode lifetimes are around 4—12 
ns [26, 27] corresponding to 7 ~ 0.01 ■ 27t GHz. 


POLARITON DECAY RATES 

Using the values above and equation 2, we plot the de¬ 
cay rates (imaginary parts of E ±) and radiative lifetimes 
of the polariton branches of an ideal Jaynes-Cummings 
system as a function of the detuning in figure 2 a and 
2b, respectively. In these plots, the red (blue) curve cor¬ 
responds to UP1 (LP1), which is more cavity-like (QD- 
like) for negative and QD-like (cavity-like) for positive 
detunings. The decay rates and thus radiative lifetimes 
exchange over a range that is determined by g and n. 

In contrast to ideal JC-systems, however, for solid- 
state systems a phonon-assisted population transfer be¬ 
tween QD excitons and cavity photons is known to exist 
[20, 28-30]. Thus, we sought to determine the true state 
lifetimes through time-resolved measurements. The fast 
cavity decay rate and strong interaction with phonons 
are expected to result in significantly shortened lifetimes. 
Hence, we performed the measurements using pulsed res¬ 
onant excitation and time-resolved detection on a streak 
camera with sub-5 ps resolution. While non-resonant ex¬ 
citation would lead to slow carrier relaxation into the 
QDs [31] our resonant scheme minimizes timing jitter in 
the excitation. To avoid exciting higher rungs of the JC 
ladder, we pumped the more QD-like polariton branch 
UP1 (as schematically illustrated by the blue arrow on 
the right side of figure la) with a pulse length that was 
carefully chosen to be shorter than the state lifetime but 
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spectrally sharp enough not to overlap with the other po- 
lariton branch. Luminescence was then observed either 
from this state directly (red arrow in figure la) or from 
the more cavity-like polariton branch after a phonon- 
assisted transfer (purple and red arrows in figure la). 
The dynamics of this process can be simulated using a 
rate equation model with four rates: radiative recombi¬ 
nation from UP1 (riypi), radiative recombination from 
LP1 (T r LP1 ), the phonon-assisted transfer rate from UP1 
to LP1 (r’f 7 '), and vice versa (r™ r ). As shown in figure 2c, 
the spectrally integrated photoluminescence intensity of 
a typical time-resolved measurement at a QD-cavity de¬ 
tuning of A = 2.8 g (including contributions from both 
decay paths) shows a mono-exponential decay, as ex¬ 
pected from the rate equation model. A fit to the data 
(red line in figure 2c) reveals a decay time of 34.9 ± 2 
ps. We note here that the spectrally integrated intensity 
decays with the same time constant as the UP1 line (see 
supplemental material for details on the rate equation 
model). 

To map out the polariton lifetimes we repeated the 
time-resolved measurement presented in figure 2c for dif¬ 
ferent detunings in the range A = 0 — llg. As shown 
in figure 2d, the resulting lifetime of UP1 increases with 
increasing detuning. However, the increase occurs much 
more slowly than indicated by the rates calculated from 
the JC model assuming no phonon interaction presented 
in figure 2b. For the detunings investigated here, LP1 is 
more cavity-like and, thus, has a much shorter radiative 
lifetime, ranging from 16 ps down to 8 ps as the detuning 
is increased. Meanwhile, the JC lifetime without phonon 
interaction for UP1 varies from 16 ps up to almost 1000 
ps, while the measured lifetime of UP1 only increases up 
to roughly 80 ps. Therefore, a strong phonon-assisted 
population transfer from UP1 to LP1 significantly short¬ 
ens the observed lifetime of UP1. 

We model this phonon-assisted population transfer us¬ 
ing an effective master equation derived in a polaron 
frame with respect to the phonon interaction [30]. This 
model fully captures the polariton dynamics and correla¬ 
tions discussed here without the need for non-Markovian 
damping. Additionally, the model is valid for our sys¬ 
tem over a wide range of temperatures, even extending 
to near absolute zero (see supplemental material for de¬ 
tails on the model and its validity). A fit to the data 
including this model is presented as a red line in figure 
2 d and produces very good overall agreement with all 
of the measured data. Since the detuning between cav¬ 
ity and QD was controlled by the lattice temperature, 
the relative temperature of each detuning was taken into 
account for these simulations. From the fit we can ex¬ 
tract the phonon-assisted population transfer rates r" r 
and r" r . They are presented in figure 2e as red and 
blue lines, respectively. Near resonance the two rates are 
very close due to the elevated temperature of T = 31 K. 
With increasing detuning both rates decrease as the po- 
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FIG. 2. (a) Decay rates and (b) resulting lifetimes of the 

polariton branches as functions of the QD-cavity detuning, 
obtained from equation 2. (c) Time-resolved spectrally inte¬ 
grated luminescence intensity at a detuning of A = 2.8 g for 
excitation of UP1. The lifetime is extracted from a mono¬ 
exponential fit (shown in red), (d) Measured and simulated 
lifetime of UP1 as a function of the detuning, (e) Phonon- 
assisted state transfer rates extracted from the lifetime at 
each detuning, (f) Spectrally and temporally resolved lumi¬ 
nescence at a QD-cavity detuning of A = 7 g. The delayed 
onset of emission from LP1 for resonant excitation of UP1 is 
consistent with phonon-assisted population transfer. 


lariton states evolve towards the bare QD and cavity 
states, where the phonon interaction cannot couple po- 
laritons through the cavity drive. The lower values of T” r 
compared to T^ r result from the different probabilities of 
phonon emission and absorption at low temperatures. 

In order to further corroborate this model of phonon- 
assisted population transfer we performed measurements 
with simultaneous spectral and temporal resolution. The 
result of a typical measurement recorded at A = 7 g is 
presented in figure 2f. Two distinct luminescence peaks 
are visible. Most importantly the onset of luminescence 
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from LP1 is delayed with respect to the onset of lumi¬ 
nescence from UP1, as population must be transferred 
before photons can be emitted. At this detuning, the ra¬ 
diative recombination rate Y r UPl is small compared to the 
strong phonon-assisted population transfer and subse¬ 
quent radiative recombination of LP1 (T^p-jJ. Hence, the 
strongest luminescence intensity is observed from LP1 
even though UP1 is resonantly excited (for details see 
supplemental material). 


PULSE-LENGTH DEPENDENT PHOTON 
BLOCKADE 

Having mapped out the detuning-dependent polariton 
lifetimes of the strongly coupled QD-cavity system we 
can now apply this knowledge to applications where the 
lifetime is critical. One of the most remarkable appli¬ 
cations of strongly coupled systems is photon blockade 
[1, 32], in which a laser pulse is tuned to the first rung 
of the JC-ladder and (due to the anharmonicity) is not 
resonant with higher transitions up the ladder. The JC 
system then transmits the portion of the incident pulse 
that couples to the system, which we observe in our cross- 
polarized reflectivity measurements. Therefore, the prob¬ 
ability of coupling a single photon through the system 
is enhanced over that of multi-photon states, which re¬ 
sults in a transmitted light beam that has a strong non- 
classical character. The fidelity of this process for the 
generation of single photons is inherently limited due to 
the linewidth of the polariton branches. Nevertheless, 
it was recently shown that in a strongly coupled system 
that is detuned by a few g 7 not only the purity but also 
the efficiency of single photon generation increases [22]. 
As we will see, photon blockade depends crucially on the 
polariton lifetimes, especially in the detuned case. 

To quantify the non-classical character of the photo¬ 
count distribution, we use its measured degree of second- 
order coherence for zero time delay [33, 34] 

s <2, (0) = < 3 > 

(n) 

computed from expectations of the measured photocount 
distribution, where n signifies a number of detections. 
Only non-classical light sources may have a second-order 
coherence g^ 2 \0 ) < 1, and g^ 2) (0) = 0 is measured exclu¬ 
sively for pulses with single-photon-like character. Owing 
to a cavity lifetime much shorter than the timing jitter of 
the single photon counters, measurements of g^(0) on 
our system can only be performed under pulsed excita¬ 
tion; this is in any case the experimental configuration 
required for on-demand applications. Hence, the choice 
of pulse length forces a compromise between frequency 
resolution (reducing the overlap of different rungs) and 
the likelihood of re-excitation of the system. In other 
words, if the laser pulses are too long the system will 
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FIG. 3. (a) Measured second-order coherence </ 2) (0) as a 

function of the laser detuning for different QD-cavity detun¬ 
ings (top to bottom) and excitation pulse lengths (left to 
right). In each case, the red vertical line indicates the bare 
cavity frequency, the green vertical line indicates the bare QD 
frequency, and the black vertical lines indicate the polariton 
frequencies, (b-c) Simulated photon-blockade values of g 1 ' 2 (0) 
as a function of the laser detuning and pulse length including 
(b) and excluding (c) the effect of phonons, both at a QD- 
cavity detuning of A = 4 g. To aid in comparison with the 
experimental results in (a), the dashed lines in (b) indicate 
the simulated anti-bunching values for photon blockade at the 
experimental pulse lengths. 


be re-excited during the interaction with a single pulse, 
reducing the non-classical character of the transmitted 
light. On the other hand, pulses with a shorter dura- 
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tion are spectrally broader, resulting in a larger overlap 
with higher rungs. With increasing detuning between 
QD and cavity resonances, the lifetime of the emitter¬ 
like (cavity-like) polariton branch increases (decreases). 
Therefore, in order to obtain the strongest photon block¬ 
ade the pulse length has to be chosen according to the 
detuning-dependent polariton lifetime. 

To test this hypothesis we performed measurements of 
g( 2 \ 0) using a Hanbury Brown and Twiss (HBT) type 
experiment [35] for different QD-cavity detunings (and 
thus polariton lifetimes) and pulse lengths. Figure 3a 
shows ( 0 ) as a function of the laser detuning for pulse 
lengths [36] of 34 ps on the left and 110 ps on the right, 
and detunings of A = 0g, A = 2g and A = 4 g from 
top to bottom. At A = 0 both pulse lengths result in a 
symmetric curve with photon-blockade dips of g ^ ( 0 ) < 1 
for laser detunings of 1 — 2 g. We measure a minimum 
g( 2 )(0) = 0.88 ±0.03 for 34 ps pulses and g^ 2 \ 0) = 0.93 ± 
0.03 for 110 ps pulses. Due to the short polariton lifetime 
of 16 ps at A = 0 the longer pulses result in an increased 
likelihood of re-excitation during the presence of a single 
pulse, which leads to a higher value of g( 2 )( 0 ) in photon 
blockade. We note that since the radiative recombination 
rates are very fast when the QD and the cavity are in 
resonance, phonons have only a minor effect in this case 
(see supplemental material for simulations and details). 

With increasing detuning the fidelity of photon block¬ 
ade increases as expected [22], and for A = 4 g (bottom 
plots of figure 3a) we measure values as low as g^ 2 \ 0) = 
0.44 ± 0.07 for 34 ps pulses and g O (0) = 0.34 ± 0.07 for 
110 ps pulses. In particular, due to the longer polariton 
lifetime at A = 4g, we observe a lower value of (/^(O) 
for the longer pulses. Here, re-excitation limits the min¬ 
imum of the second-order coherence for the long pulses, 
while spectral overlap with transitions to higher rungs 
limits the minimum for the short pulses. 

To explain this finding we performed quantum opti¬ 
cal simulations using the Quantum Optics Toolbox in 
Python (QuTiP) [37] based on the Quantum Regres¬ 
sion Theorem [22], The results of these simulations are 
presented in figure 3b-c. The figures show </ 2 )(0) as 
a function of the laser detuning and pulse length for 
A = 4 g, both including (figure 3b) and excluding (figure 
3c) phonon-assisted population transfer (see supplemen¬ 
tal material for details on the simulations and additional 
representations). In both cases, with increasing pulse 
length the photon-blockade dip narrows down (as seen in 
the experiment) due to the improved spectral resolution 
of longer pulses. Moreover, the minimum value of (/ 2 f (0) 
decreases due to the decreased overlap with higher rungs 
(as well as the other polaritonic branch) before increas¬ 
ing due to re-excitation for too-long pulses. As can be 
seen by comparing figure 3b with figure 3c, the phonon- 
induced reduction of the polariton lifetime means that 
the deepest photon-blockade dip is observed at a shorter 
pulse length when phonons are included, and the corre¬ 


sponding minimum g^(0) values are a little higher. The 
simulations are in very good overall agreement with the 
measurements presented in figure 3a, and they highlight 
the importance of taking plronon-assisted transfer into 
account when selecting the optimal pulse length. 

For an excitation between the polariton branches, a 
peak with <je 2 )(0) > 1 is observed for all detunings and 
pulse lengths. This phenomenon is known as photon- 
induced tunneling [1], where anharmonicity of the JC 
ladder results in bunching. However, the bunching val¬ 
ues peak at the position where the probability of zero 
detections is highest and simulations of the plrotodetec- 
tion distributions barely deviate from those of a coherent 
state with a similar count rate [22], Therefore, in this re¬ 
port we focus on the antibunching regime. 

COHERENTLY-EXCITED SINGLE-PHOTON 
SOURCE 

While we have seen above that for detuned photon 
blockade the presence of exciton-phonon coupling slightly 
decreases the fidelity, this coupling also allows us to inves¬ 
tigate more sophisticated schemes for on-demand single¬ 
photon generation. In particular, the population of UP1 
can be coherently controlled if the pulse length is cho¬ 
sen shorter than the state lifetime but spectrally narrow 
enough to avoid exciting higher transitions up the JC lad¬ 
der. Due to coherent scattering of the laser and imperfect 
suppression of the laser reflection from the sample sur¬ 
face, it is difficult to observe Rabi oscillations directly. 
However, phonon-assisted emission from LP1 following 
excitation of UP1 occurs at a different frequency. In ad¬ 
dition, for detunings in the range of A = 5 — 10 <?, T r UP1 
is strongly Purcell-suppressed while T£ P1 is very fast. 
Most strikingly, Vy is proportional to g 2 [30] and, thus, 
very efficient for a strongly coupled system. Therefore, 
even when resonantly exciting UP1 for these detunings 
most luminescence is emitted from LP1. The result of 
such an experiment is presented in figure 4a, which shows 
the emission intensity from LP1 for a resonant excitation 
of UP1 with 16 ps pulses at a detuning of A = 4.8 g. 
Clear Rabi oscillations are observed and fitted (red line 
in figure 4a) using a damped oscillation superimposed 
with a very weak background that is linear with the ex¬ 
citation power, due to a finite leakage of the excitation 
laser into the detection channel. This power-dependent 
damping is a direct consequence of the polariton-phonon 
interaction. With increasing laser power the polariton 
levels hybridize, mixing higher rungs with UP1. This 
increases the number of available damping channels via 
the polariton-phonon coupling and results in a damped 
oscillation. 

To investigate single photon generation in this config¬ 
uration we measure t/( 2 )(0) while exciting UP1 with a tt 
pulse and detecting LP1 emission. The result of this ex- 



6 



FIG. 4. (a) Rabi oscillation observed upon exciting UP1 and 
detecting LP1. (b) Measurement of the second-order coher¬ 
ence for exciting UP1 with a n pulse and detecting LP1 yield¬ 
ing fl (2) (0) = 0.03 ± 0.01 


periment is presented in figure 4b and shows almost per¬ 
fect antibunching. By integrating the area of the peaks 
we obtain a value of g^ 2 \ 0) = 0.04 ± 0.01. Fitting the 
data with a series of Gaussian peaks is a little less sen¬ 
sitive to noise from the dark counts of the detectors and 
yields g^{ 0) = 0.03±0.01 (see supplemental material for 
the fit and a closer examination of the data). These re¬ 
sults rival the best values obtained from QDs, while gain¬ 
ing advantages provided by the optical cavity (e.g., the 
very short lifetime that enables a high photon generation 
rate and the potential for integration into on-chip optical 
networks). Hence, our findings clearly demonstrate the 
potential of exploiting the efficient exciton-phonon cou¬ 
pling of strongly coupled QD-cavity systems for single 
photon generation. 


SUMMARY AND CONCLUSIONS 

The presence of phonons represents one of the key dis¬ 
tinguishing features of the solid state. Our investigations 
have laid the groundwork for understanding their role 
in the dynamics of strongly coupled systems. We have 
mapped out the detuning-dependent polariton lifetime of 
a strongly coupled QD-cavity system and extracted the 
spectrum of the polariton-to-phonon coupling for QD- 
cavity detunings up to 11 g. We have shown that in or¬ 
der to obtain high-fidelity photon blockade a pulse length 
must be chosen that depends on the QD-cavity detun¬ 
ing. Finally, we have achieved direct coherent control of 
polariton states and provided an unprecedented demon¬ 
stration of how the efficient coupling to phonons can be 
exploited for high-fidelity single photon generation. Our 
findings suggest that the influence of the efficient cou¬ 
pling to phonons should be considered for all applications 
of strongly coupled solid-state systems. Furthermore, 
we expect that our contribution towards understanding 
strongly coupled systems in the solid state will play a vi¬ 
tal role in utilizing QD-cavity platforms for novel physics 


such as higher-order non-classical light generation [3]. 
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SUPPLEMENTAL MATERIAL 


METHODS 

Sample fabrication 

The MBE grown structure consists of a ~ 900 nm thick 
Al 0 .sGa 0 . 2 As sacrificial layer followed by a 145 nm thick 
GaAs layer that contains a single layer of InAs QDs. 
Our growth conditions result in a typical QD density 
of 60 — 80 /mi -2 . The photonic crystals were fabricated 
using 100 keV e-beam lithography with ZEP resist, 
followed by reactive ion etching and HF removal of the 
sacrificial layer. The photonic crystal lattice constant 
was a = 246 nm and the hole radius r ss 60 nm. The 
cavity fabricated is a linear three-hole defect (L3) cavity. 
To improve the cavity quality factor, holes adjacent to 
the cavity were shifted [39, 40]. 

Optical spectroscopy 

All optical measurements were performed with a liquid 
helium flow cryostat at temperatures in the range 
10 — 40 K. For excitation and detection a microscope 
objective with a numeric aperture of NA = 0.75 was 
used. Cross-polarised measurements were performed 
using a polarising beam splitter. To further enhance the 
extinction ratio, additional thin film linear polarisers 
were placed in the excitation/detection pathways and a 
single mode fibre was used to spatially filter the detec¬ 
tion signal. Furthermore, two waveplates were placed 
between the beamsplitter and microscope objective: 
a half-wave plate to rotate the polarisation relative 
to the cavity and a quarter-wave plate to correct for 
birefringence of the optics and sample itself. 

Autocorrelation measurements 

Second-order autocorrelation measurements were per¬ 
formed using a Hanbury Brown and Twiss (HBT) setup 
consisting of one fibre beamsplitter and two single 
photon avalanche diodes. The detected photons were 
correlated with a PicoHarp300 time counting module. 

Simulations 

The simulations of g 2 (0) presented in figure 3 of the 
main text were performed using the Quantum Optics 
Toolbox in Python (QuTiP) [37]. Here, we used a 
Jaynes-Cummings Hamiltonian evolution subject to 
dissipation terms of the Lindblad form. These terms 
include those due to the cavity decay and the phonon- 
induced cavity—^exciton and exciton—i-cavity processes, 
with rates derived from the phonon model discussed 
below. The simulations are time-dependent, in that the 





FIG. 5. (a) Calculated dynamics of the state population of 
UP1 and LP1 after an excitation of UP1 at A = 3 g. (b) 
Calculated dynamics of the emission intensity from UP1 and 
LP1. (c) Calculated dynamics of the spectrally integrated 
emission intensity. 


cavity is driven with a Gaussian pulse of experimentally 
extracted pulse-width. The only free parameter is the 
laser power. Finally, the calculation of the pulse-wise 
second order coherence is exact, in that it is the sum 
of all two-time photon-photon correlations of cavity 
photons. These correlations are calculated via the 
quantum regression theorem [22]. 
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RATE EQUATION MODEL FOR DESCRIBING 
THE SYSTEM DYNAMICS 


As discussed in the main text, we use a rate equa¬ 
tion model with two states and four rates to describe 
the system dynamics: radiative recombination from UP1 
(r^pi), radiative recombination from LP1 (T r LP1 ), the 
phonon-assisted transfer rate from the UP1 to LP1 (r^ 7- ), 
and vice versa (r” r ). The resulting rate equation model 
then reads: 


{Pupi{t)\ 
dt \PLPi(t)) 

( /rr _i_ T^nr\ ~pnr 

\ L UP 1 + 1 / ) 1 r 

Y'nr (j^r _±_ ~pnr\ 

1 f V 1 LP1 i 1 r ) 


( Pupi(t)\ 
\PLPl(t)J 


( 4 ) 


where Pjjpi(t) and Pipi(t) describe the populations 
of the states UP1 and LP1, respectively. From solving 
these equations we obtain the time dependence of the 
state populations. 

An example is presented in figure 5a that shows the 
populations of UP1 and LP1 at a detuning of A = 3g 
in red and blue for excitation of UP1 and using the val¬ 
ues for the four rates as obtained in the main text for 
this detuning. After a short time the populations quasi- 
thermalize and decay with an almost perfect exponential 
shape. The resulting intensity of emission from these two 
states is presented in figure 5b using the same colour cod¬ 
ing as figure 5a and shows the same behaviour. The spec¬ 
trally integrated emission intensity is presented in figure 
5c. Importantly, after the fast quasi-thermalization the 
emission intensities from the two states (figure 5b) decay 
with the same time constant. Therefore, the spectrally 
integrated intensity (figure 5c) decays with the same time 
constant. This allows extraction of the lifetime of UP1 
from temporally analysing the spectrally integrated lu¬ 
minescence intensity. In particular, the streak camera 
that was used to record the temporally resolved spectra 
has a lower spectral resolution than the spectrometer and 
CCD that was used to record the time integrated spectra 
as presented in figure 1 of the main text. Therefore, for 
detunings smaller than A « 5 g analysing the spectrally 
integrated intensity is the only option as the two time 
traces cannot be spectrally resolved. We note here that 
these finding hold for all investigated sets of parameters. 


tral function is used which is given by 

J(w) = a p u 3 exp ( 5 ) 

and describes the interaction between electrons and the 
logitudinal acoustic phonons via deformation potential 
coupling. This function uses two parameters a p and u>i 
The latter one corresponds to a high frequency cutoff due 
to the electronic localization length. The fit of the de¬ 
tuning dependent polariton lifetime presented in figure 2 
of the main text resulted in the values a p = 3 ps 2 and 
u>b = 0.22 meV. However, we note that these values vary 
quite a lot throughout literature. This can be understood 
from the fact that J(ui) is derived for a spherical quan¬ 
tum dot with similar extents for electron and hole wave- 
functions. Moreover, experiments typically cover only a 
limited frequency range and therefore can be fitted with 
a range of parameters. 

We now discuss the validity of our application of the 
effective phonon master equation. Given that the EPME 
results from a truncation of the cavity-exciton transfer 
process to first order, it is incapable of modelling the 
asymmetry in the forward and reverse rates that physi¬ 
cally arises due to spontaneous phonon emission. How¬ 
ever, at our experimental bath temperatures the phonon 
occupancies are elevated such that the EPME only in¬ 
correctly estimates these transfer rates by a few per¬ 
cent - well within experimental error. Furthermore, even 
at low temperatures the EPME is capable of accurately 
capturing the dynamics of a typical QD-cavity system. 
Here, the radiative decay rate dominates the state decay 
near zero detuning, where the asymmetry in the phonon- 
induced transfer rates is the largest. Away from zero 
detuning the EPME correctly captures the influence of 
spontaneous phonon emission because the lowest-order 
polaritonic states have a clear excitonic or photonic char¬ 
acter that can be matched to the transfer terms. There¬ 
fore, we are justified in ignoring these small inaccuracies 
and modeling our system dynamics with the EPME. We 
also note here that spectral diffusion does not affect the 
ultrafast dynamics described here as it occurs on a much 
slower timescale [12]. 


PHONON MODEL 


To model the polariton-phonon coupling we used a 
model that utilizes an effective master equation derived 
in a polaron frame with respect to the phonon interaction 
(also known as the effective phonon master equation, or 
EPME) [30]. In this model a characteristic phonon spec- 
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FIG. 7. Comparison of measured and simulated g 2 ( 0) as func¬ 
tions of the laser detuning for a QD-cavity detuning of A = 4 g 
with (a) 34 ps and (b) 110 ps long pulses. 


FIG. 6. Simulated g 2 ( 0) as a function of the laser detun¬ 
ing and pulse length for a QD-cavity detuning of A = 0 (a) 
including and (b) excluding phonons 

PHOTON BLOCKADE SIMULATIONS 

In the main text we presented simulations of <? 2 (0) 
as a function of the pulse length and laser detuning 
for A = 4 g with and without phonons. As mentioned 
there, we also performed similar simulations for A = 0 g. 
The result of these simulations is presented in figure 6a 
with phonons and figure 6b without phonons. In both 
cases the best photon-blockade is observed for relatively 
short pulses around 40 ps which can be understood from 
the short polariton lifetime at resonance which results 
in strong re-excitation of the system for longer pulses. 
Moreover, compared to the situation at A = 4 g (main 
text figure 3b and 3c) the difference in the optimum pulse 
length as well as observed value of g 2 {0 ) for the cases 
with and without phonons is much smaller. This results 
from the fact that the radiative decay rate of the po¬ 
lariton branches at resonance is much faster than the 
phonon mediated population transfer. Thus, the effect 
of phonons is much less pronounced. Regardless, since 
photon-blockade has a much higher efficiency at a detun¬ 
ing of a few g compared to the resonant case, the impact 
of phonons on the polariton lifetime has to be considered 
carefully in order to optimize the experimental conditions 
for single photon generation by photon blockade. 


We also present a direct comparison between the mea¬ 
surements and simulations of g^ (0) from figure 3 of the 
main text. To this end, we replotted in figure 7 the 
measurements for A = 4 g (from figure 3a of the main 
text) adjacent to simulations (from figure 3b of the main 
text) using the same representation. While there is some 
jitter in the x-values of the measurements due to the 
limited resolution of our frequency filter and tempera¬ 
ture drift during the measurements, the simulations re¬ 
produce the experiments very well. In particular, the 
photon-blockade dip is quantitatively well reproduced. 
The simulations also reproduce the experimental trend 
that the bunching peak is higher for longer pulses. How¬ 
ever, the bunching values should be loosely interpreted 
given that they occur at the position where the probabil¬ 
ity of zero detections is highest (as discussed in the main 
text). Thus, the bunching values are strongly altered by 
effects such as blinking of the QD or an imperfect sup¬ 
pression of the excitation laser. 

PHONON-ASSISTED SINGLE PHOTON 
GENERATION 

In the main text, we presented the measurement of 
g! 2 ^(0) for phonon-assisted single photon generation on 
a linear scale. Here, we represent the data on a semi- 
logarithmic scale in figure 8a to further visualise the sin¬ 
gle photon purity. A fit to the data with a series of peaks 
is presented as a red line and a magnification of the cen- 
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ter peak is presented in figure 8b. Due to the short po- 
lariton lifetime compared to the repetition rate, a flat 
background resulting only from dark counts can be seen. 
This can be used to fit the average number of dark count 
coincidences, which due to the high temporal resolution 
of our measurement is only 1.88. 


FIG. 8. (a) Semi-logarithmic representation of the phonon- 
assisted single photon generation presented in figure 4b of the 
main text. A fit with a series of peaks is presented as a red 
line (b) Zoom on the center peak. 







































